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We briefly sketch the methods for a numerically stable evaluation of tensor one-loop integrals that have been 
used in the calculation of the complete electroweak one-loop corrections to e + e~ — > 4fermions. In particular, 
the improvement of the new methods over the conventional Passarino-Veltman reduction is illustrated for some 
4-point integrals in the delicate limits of small Gram (and other kinematical) determinants. 
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1. INTRODUCTION 

At the LHC and ILC, many interesting pro- 
cesses involve more than four external particles. 
Such many-particle reactions often proceed via 
resonances that subsequently decay, or they rep- 
resent a background to resonance processes. A 
thorough description of such processes requires 
the evaluation of strong and electroweak radiative 
corrections at least in next-to-leading order. The 
most complicated part in such calculations con- 
cerns the numerically stable evaluation of the one- 
loop tensor integrals of the virtual corrections. 

For processes with up to four external particles 
the classical Passarino-Veltman (PV) reduction 
[1] , which recursively reduces tensor to scalar in- 
tegrals, is sufficient in practically all cases. This 
scheme, however, involves Gram determinants in 
the denominator, which spoil the numerical sta- 
bility if they become small. With up to four ex- 
ternal particles this happens only near the edge 
of phase space (forward scattering, thresholds). 
With more than four external particles, Gram de- 
terminants also vanish within phase space, and 
methods are needed where Gram determinants 
can be small but still non-zero. Several solutions 
to this problem have been proposed in recent 
years, but not many of them have proven their 
performance in complicated applications yet. For 
references and descriptions of some methods al- 
ternative to ours, we refer to Refs. [2,3]. 

In this article we briefly describe the meth- 
ods used in the calculation of the complete elec- 



troweak 0(a) corrections to 



4 fermions 



[4], which constitutes the first established one- 
loop result for a 2 — > 4 particle reaction. All rele- 
vant formulas can be found in Refs. [2,5], here we 
only sketch their structure. Moreover, we illus- 
trate the improvement of the new methods over 
the conventional PV reduction for some 4-point 
integrals in delicate kinematical limits. 



2. THE GENERAL CONCEPT 

One-loop tensor integrals can be naturally 
grouped into three categories, which we have 
treated in completely different ways: 

(i) For 1- and 2-point integrals of arbitrary ten- 
sor rank, numerically stable analytical expres- 
sions are presented in Ref. [2] (see also Ref. [1]). 

(ii) For 3- and 4-point tensor integrals. PV re- 
duction [1] is applied for "regular" phase-space 
points where Gram determinants are not too 
small. For the remaining problematic cases spe- 
cial reduction techniques have been developed [2] . 

One of the techniques replaces the standard 
scalar integral by a specific tensor coefficient that 
can be safely evaluated numerically and reduces 
the remaining tensor coefficients as well as the 
standard scalar integral to the new basis inte- 
grals. In this scheme no dangerous inverse Gram 
determinants occur, but inverse modified Caylcy 
determinants instead. The procedure is related to 
the fully numerical method described in Ref. [6]. 

In a second class of techniques, the tensor co- 
efficients are iteratively deduced up to terms that 
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are systematically suppressed by small Gram or 
other kincmatical determinants in specific kinc- 
matical configurations. The numerical accuracy 
can be systematically improved upon including 
higher tensor ranks. A similar idea, where ten- 
sor coefficients are iteratively determined from 
higher-rank tensors has been described in Ref. [7] 
for the massless case. 

These methods are sketched in more detail and 
illustrated in 4-point examples below. 

(iii) For 5- and 6-point integrals, direct reduc- 
tions to 5- and 4-point integrals, respectively, 
are possible owing to the four-dimensionality of 
space-time. For scalar integrals such a reduc- 
tion was already derived in the 1960s [10]. In 
Rcfs. [2,5] we follow basically the same strategy to 
reduce tensor integrals, which has the advantage 
that no inverse Gram determinants appear in the 
reduction. Instead modified Cayley determinants 
occur in the denominator, but we did not find nu- 
merical problems with these factors. A reduction 
similar to ours has been proposed in Ref. [8] . 

We would like to stress two important features 
of our approach. 

(i) The methods are valid for massive and mass- 
less cases. The formulas given in Refs. [2,5] are 
valid without modifications if IR divergences are 
regularized with mass parameters or dimension- 
ally. 1 Finite masses can be either real or complex. 

(ii) The in/out structure of the methods is the 
same as for conventional PV reduction, i.e. no 
specific algebraic manipulations are needed in ap- 
plications. Therefore, the whole method can be 
(and in fact is) organized as a numerical library 
for scalar integrals and tensor coefficients. 

We conclude this overview with some com- 
ments resulting from our experience collected in 
the treatment of a full 2^4 scattering reaction. 

(i) For a specific point in a multi-particle 
(multi-parameter) phase space it is highly non- 
trivial to figure out which of the various methods 
is the most precise. It seems hopeless to split the 
phase space into regions that are dedicated to a 
given method. Therefore, we estimate the accu- 
racy for the different methods at each phase-space 



point and take the variant promising the highest 
precision. In one approach, we estimate the num- 
ber of valid digits based on the expected accuracy 
of the expansions and possible numerical cancel- 
lations before the evaluation of the coefficients. 
In a second approach, the accuracy of the PV 
method is estimated by checking symmetries and 
sum rules for the coefficients. If the estimated 
precision is not satisfactory, the seminumerical 
method is taken; if this is still not satisfactory, 
an expansion is used. 

(ii) In a complicated phase space it may hap- 
pen that none of the various methods is perfect 
or good in some exceptional situations. Usually 
the corresponding events do not significantly con- 
tribute to cross sections. This issue can only be 
fathomed in actual applications. To be on the safe 
side, we employ the two independent "rescue sys- 
tems" with different advantages and limitations. 

(iii) In view of this, figures as shown below are 
nice illustrations, but should always be taken with 
a grain of salt. No matter how many of such 
figures are shown, they will never be exhaustive, 
so that no quantitative conclusions on the overall 
precision of methods can be drawn. 

In summary, the performance (speed, stability, 
etc.) of any method can only be estimated and 
proven in non-trivial and realistic applications. 

3. 4-POINT INTEGRALS FOR SMALL 
GRAM DETERMINANTS 

In the following we sketch the methods of 
Ref. [2] for the 4-point tensor integrals of rank P, 
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N k = (q + pkf - m\ + ic, po = 0, (1) 
which are decomposed into covariants as follows, 



1 For the method of Ref. [5] , this has been shown in Ref. 
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E fe"X + S"X + s"X ) Ax* , (2) 
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and so on for higher rank. The conventional PV 
reduction [1] recursively expresses the tensor co- 



efficients Di 



in terms of the scalar 4-point 



integral Dq and 3-point integrals of lower rank. 
To this end, the defining tensor integrals are con- 
tracted with an external momentum p^ 1 or with 
the metric tensor g^ 1 ^ 2 , leading to the factors pkq 
or q 2 in the numerator, respectively. Rewriting 
these factors according to 

2p k q = N k -N - f k , f k =Pk~ m 2 k + ml, 
q 2 = N Q + m 2 , (3) 

the N terms can be used to cancel propagator de- 
nominators. This procedure yields the basic PV 
relations that connect contracted 4-point tensor 
integrals of rank P with lower-rank 4- and 3-point 
integrals. The PV solution expresses the rank P 
4-point coefficients in terms of lower-rank 4- and 
3-point coefficients. In each step P — > (P— 1) the 
inverse of the Gram matrix 

2pipi 2pxp 2 2pxp 3 \ 
2p 2 pi 2p 2 p2 2p 2 p 3 (4) 
2p 3 pi 2p 3 p 2 2p 3 p 3 J 

occurs, which causes the above-mentioned nu- 
merical problems if the determinant \Z\ becomes 
small. In our alternative methods also the matrix 



/ 2m\ 
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fi h h 
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(5) 



and its inverse appear. The vanishing of the mod- 
ified Cayley determinant \X\ corresponds to nec- 
essary conditions for true (Landau) singularities 
in a Feynman diagram. 

We adopt the conventions that indices i with 
a hat should be omitted and that and X^ 
denote the minors (i.e. determinants of subma- 
trices where row i and column j are discarded) 
of the matrices Z and X, respectively. Similarly 
Zij,kt denotes generalized minors where two rows 
and two columns are discarded in Z. For the 
precise definitions, full formulas, and detailed de- 
scriptions of algorithms we refer to Ref. [2] . Here 
we only highlight the general structure. 



3.1. Seminumerical method for small \Z\ 

The PV relations can be arranged as follows, 
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3-point integrals, (6) 



Er=2 S 2ir D 
\J2r=2 6 ^r D 



which expresses tensor coefficients of rank P and 
(P — 1) in terms of coefficients of rank P and 
(P + 1) that involve one metric tensor more in 
the corresponding covariants. Thus, these rela- 
tions can be used to recursively express all 4-point 
coefficients up to rank P (including the scalar in- 
tegral) in terms of the coefficient -D0...0 of rank 
2P and 3-point integrals. For rank 4 or higher 
the coefficient -D0...0 has the property that its 
Feynman-parameter integral involves a logarith- 
mic integrand without any denominator and is, 
thus, well suited for numerical integration. 

Figure 1 illustrates the reliability of the seminu- 
merical method for an example where the Gram 
determinant \Z\ becomes small (x — > 0) and 
therefore the PV algorithm breaks down. Away 
from the tricky region, i.e. for increasing x (and 
thus \Z\), there is a transition region in which PV 
reduction and the seminumerical method agree 
within good accuracy until the latter becomes 
less precise than the PV method. This is due to 
the error of the numerically evaluated coefficient 
-D0...0 which enters with a prefactor \Z\, and is 
thus suppressed for small \Z\. 

The seminumerical method is limited to cases 
in which \X\ is not too small. If necessary, the 
case of both \Z\ and \X\ small, where PV reduc- 
tion does not work either, can be cured with an 
expansion described below. For 3-point integrals, 
the soft or collincar singular cases have the prop- 
erty that \X\ =0 exactly, so that the seminu- 
merical method is not applicable. However, these 
special cases are simple enough to be dealt with 
analytically as described in App. B of Ref. [2]. 

3.2. Expansion for small \Z\ 

The PV relations can be rewritten as 
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Full diagram: 
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Figure 1. A typical example for 4-point integrals with small \Z\ (x — > 0). The full diagram and the relevant 
subdiagram are given above; absolute predictions (in arbitrary units) for some tensor coefficients, relative 
deviations from PV reduction, and the kinematic specifications are shown below. 
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Figure 2. An example for 4-point integrals with both \Z\ and \X\ small (x — ► 0). Details as in Fig. 1. 
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+ 2m z Z k iD ll ... lp - \Z\D kHl ... ip 

+ 3-point integrals, (7) 

which express the coefficients D^.ip and 
DQ0h...i P hi terms of lower-rank 4- and 3-point 
coefficients and higher-rank 4-point coefficients 
dressed with a prefactor \Z\. For small \Z\, these 
relations can, thus, be used for an iterative deter- 
mination of all 4-point coefficients from 3-point 
coefficients up to terms that are suppressed by 
some powers of \Z\. For a specific coefficient the 
iteration can be made arbitrarily precise upon in- 
cluding higher and higher ranks in the iteration. 

Figure 1 demonstrates the reliability of the 
method for the same configuration as for the 
scminumcrical method of the previous section. 
For increasing \Z\ (and x) the expansion becomes 
less and less precise, because the missing terms 
that are suppressed by some power of \Z\ grow. 
Again there is a region where all methods yield 
decent results. 

The discussed expansion is limited to the case 
where Xqj and Z^i are not too small for at least 
one set of indices j, k, I. If all Xqj are small, then 
\X\ is small, too. This case is considered in the 
next subsection. The case in which all Zki axe 
small is further elaborated in Ref. [2]. 

3.3. Expansion for small \Z\ and |X| 

The PV relations can again be rewritten as 

p 

y^fcy-PpQ^ ,..i r ...ip x Xs,„D, ...... - \Z\D kil ... lp 

r=l 

+ 3-point integrals, 

1),, = const. X ZijD? l Qi 1 ...i p 
P N-l 

_2 EE %«)(jv)/™- D ooii...i r ...; P 



+ XQjDu 1 ...i p + 3-point integrals, (8) 

which expresses the coefficients D^ . ip and 
Dooii...ip in terms of lower-rank 4- and 3-point 
coefficients and higher-rank 4-point coefficients 
dressed with a prefactor \Z\ or Xoj. If \Z\ and 
all Xoj are small, these relations can be used to 
itcratively determine all 4-point from 3-point co- 
efficients up to suppressed terms. Compared to 
the case of the previous subsection, where only 
\Z\ was considered small, the tensor rank grows 
faster with each iteration. 

Figure 2 shows the reliability of the method for 
a specific example, revealing again the breakdown 
of PV reduction for sufficiently small x (and thus 
\Z\ and \X\), while the expansion becomes less 
precise for increasing x. 

The expansion method fails if either all Zki r 
or all Xij are small. Possible treatments of these 
exceptional cases are also described in Ref. [2]. 
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